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ABSTRACT 


Fourier  spectrometry  is  at  a  disadvantage  compared  to 
conventional  scanners  insofar  as  larger  dynamic  range  is  re¬ 
quired  for  recording  and/or  transmitting  the  signal.  Various 
procedures  for  minimizing  the  disadvantage  are  evaluated.  A 
simulated  experimental  comparison  is  made  between  signal 
compression  techniques  and  chirping  (distribution  of  the  central 
fringes  of  an  interferogram) .  The  comparison  yields  approx¬ 
imately  a  15  dB  improvement  for  both  procedures,  though  each 
leave' somewhat  different  residual  distortions  of  the  spectrum. 
The  emphasis  is  given  to  the  chirping  procedure,  since  this 
is  less  well  known,  and  ascendence  over  its  problems  provides 
greater  insight  into  Fourier  spectrometry. 


SECTION  I 


INTRODUCTION 


A  recurring  problem  in  the  implementation  of  Fourier 
spectrometry  has  been  the  high  precision  required  for  the 
recording  of  interferograms.  This  results  from  the  fact 
that  interferograms  usually  have  a  much  greater  dynamic 
range  (peak  signal-to-noise  ratio)  than  the  spectra  to 
which  they  correspond.  Although  the  problem  is  most  severe 
in  systems  where  telemetry  xs  employed,  dynamic  range  can 
pose  a  limitation  in  laboratory  and  field  experiments  as 
well. 
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SECTION  II 

DYNAMIC  RANGi  REQUIREMENTS 


Any  device  used  to  transmit,  record  or  in  any  way 
process  an  '  ' terferogram  must  have  a  dynamic  range  capabil¬ 
ity  at  least  equal  to  the  dynamic  range  of  the  inter ferogr am 
if  no  information  is  to  be  lost.  The  expression1  relating 
interferogram  dynamic  range  to  the  mean  spectral  signal - 
to-noise  ratio  is 

Dynamic  range  =  I  /M  =  -JF  Jav  /Av  ’  ~  m 

°  I  *  *  s  m  M 

s 

where  Iq  is  the  peak-to-peak  amplitude  of  the  envelope  of 

the  central  fringes  of  the  interferogram,  S  the  mean  spectral 

amplitude,  N  che  number  of  spectral  resolution  elements,  and 

M];  and  Mg  the  root  mean  square  noise  on  the  interferogram 

and  spectrum  respectively.  Avg  and  Avm  are  respectively 

the  frequency  range  of  the  spectral  information  and  that  of 

the  noise.  Since  the  system  bandwidth  need  not  be  greater 

than  that  necessary  to  process  the  signal,  Av  /Av  will  be 

s  m 

considered  unity.  Figure  1  depicts  equation  (1)  for  a  range 
of  resolving  powers. 

Typical  dynamic  range  requirements  have  ranges  from  about 
3  4 

10  to  10  .  The  spectrometer  described  by  the  Connes2  resolves 
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on  the  order  of  10  'spectral  elements  and  can  record  inter- 

0 

ferograms  with  a  dynamic  range  o^  10  .  From  figure  1  we 
conclude  that  the  mean  signal-to-noise  ratio  possible  with 
that  instrument  is  well  in  excess  of  1000. 

In  cases  where  the  dynamic  range  is  limited  by  digitiz¬ 
ing  noise  equacion  (1)  may  not  be  applicable.  The  reason  is 
that  digitizing  noise  can  affect  interferograms  in  a  highly 
systematic  manner.  For  example,  in  a  system  limited  by 
digitizing  noise  it  is  always  possible  that  some  fringes  lie 
below  the  digitizing  threshold.  The  consequences  can  range 
from  premature  truncation  oi  the  interferogram  to  the  pro¬ 
duction  of  blank  segments.  Premature  truncation  results  in 
a  loss  of  spectral  resolving  power  while  the  presence  of 
blank  segments  can  totally  destroy  the  instrumental  profile. 

the  other  hand,  if  the  weak  fringes  happen  to  lie  just  at 
the  digitizing  threshold,  then  they  are  unduly  amplified. 

To  prevent  such  difficulties  it  is  essential  that  the  random 
inter ferogram  noise  be  greater  than  the  least  significant 
digit. 

Equation  (1)  determines  the  dynamic  range  necessary 
where  a  single  ordinary  interferogram  is  used  to  compute  the 
spectrum.  Interferogram  processing  devices  need  not  process 
only  one  interferogram,  per  spectrum,  and  interferograms  can 
be  of  types  other  than  ordinary.  The  following  sections  deal 
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with  these  topics. 
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SECTION  III 


BANDWIDTH  MANIPULATION 


Dynamic  Range  -  Bandwidth  Trade-off: 

A  situation  is  often  encountered  in  which  an  interfero- 
gram  processing  device  of  limited  dynamic  range  posseses  a 
surplus  of  available  bandwidth.  In  such  cases  full  use  of 
the  available  bandwidth  will  substantially  reduce  dynamic 
range  requirements.  Before  presenting  this  procedure,  a 
discussion  of  interf erogram  bandwidth  (including  a  brief  re¬ 
view  of  Fourier  spectrometer)  is  in  order 

Interf erogram  Bandwidth; 

Figure  2  represents  a  typical  Fourier  spectrometer 

k 

configuration.  Assume  that  the  moving  mirror  is  uniformly 
translated  such  that  in  time  T  the  path  difference  between 
the  two  interferometer  beams  has  varied  around  zero  over  a 
range  +_  At/2.  For  each  spectral  component  of  the  incident 
radiation  of  frequency  v  cm  \  the  detector  experiences 
vAt  cycles  of  cosinusoidal  modulation  in  time  T.  An  electri¬ 
cal  frequency  f  corresponding  to  v  cm  ^  is  thus  produced  by 
the  detector  according  to 

f  =  vAt/T  Hz  .  (2) 

* 

The  difference  between  continuous  scanning  and  step-by-step 
scanning  is  slight  in  this  context.  Refer  to  reference  1 
for  a  discussion  of  this  topic. 
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F«XE.D  KURROR 


Fi$.  ’  Mickeleon  interferometer  fer  Fourier  fpeetrometrj. 
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Since  the  spectral  resolution  is  given  by 


5v  =  ~~  cm 

At 


the  bandwidth  of  the  inter ferogr am  is 


(3) 


Af 


v 

6vT 


(4) 


where  N  is  the  number  of  spectrum,  resolution  elements. 
Required  Measurement  Time: 

If  the  interferogram  noise  n^  in  equation  (1)  is  white 
detector  noise  it  can  be  determined  from  the  noise  equivalent 
power  of  the  detector  by 

n  =  N.E.P.  (Af ) 2  watts  (5) 


where  the  noise  bandwidth  is  limited  to  that  of  the  signal. 
The  dynamic  range  therefore  becomes 


I  /N.E.P.  (Af) 
o 


By  substituting  into  equation  (1),  the  signal-to-noise 
ratio  obtainable  in  measurement  time  T  is  given  by 

-  I  J~T~' 

S  _  o  V  _ 

n  “  J.E.P.  N 

s 
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Multiple  Scanning: 

If  the  entire  measurement  time  required  for  a  satis¬ 
factory  signal- to-noise  ratio  is  alloted  to  a  single  mter- 
ferogram,  Af  is  minimum,  and  maximum  dynamic  range  is 
required.  The  measurement  time  may  however  be  used  to  observe 
k  inter ferograms  each  of  duration  T/k.  For  this  situation 
Af  becomes  kN/T  and  the  dynamic  range  is  reduced  by^k 
The  data  reduction  procedure  requires  the  averaging  of  k 
interferograms  in  order  to  obtain  the  required  signal- to- 
noise  r?tio.  Although  the  composite  inter ferogram  will  have 
the  original  large  dynamic  range,  it  need  only  reside  in  a 
computer  memory  where  dynamic  range  is  readily  available. 

Sufficient  bandwidth  must  be  freely  available  for  this 
procedure  to  be  appropriate,  since  it  does  not  pretend  to  use 
the  bandwidth  most  efficiently.  Thus,  it  only  improves  the 
performance  of  '\nalog  systems,  since  dynamic  range  and  band¬ 
width  are  inseparable  on  digital  channels. 

Examples : 

Multiple  scanning  can  be  used  to  reduce  dynamic  range 
where  the  limitation  is  imposed  by  analog  telemetry  of  tape 
recording.  These  devices  are  typically  limited  to  a  dynam¬ 
ic  range  of  about  100:1  and  are  capable  of  operating  over 
large  bandwidths. 
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As  a  typical  case  consider  a  Fourier  spectrometer  re¬ 
solving  1000  spectral  elements.  Suppose  that  a  signal-to- 
noise  ratio  of  100  is  desired  and  that  the  data  are  to  be 
recorded  on  an  analog  tape  recorder  with  a  dynamic  range  of 
100.  Suppose  further  that  equation  (6)  indicates  an  observing 
time  of  1000  seconds  is  required  for  the  particular  source 
and  detector.  Figure  1  indicates  that  a  dynamic  range  of 
about  3000  is  required  and  equation  (4)  indicates  a  band¬ 
width  of  1  Hz.  To  properly  record  the  signals,  the  observ- 

2 

ing  time  should  be  divided  among  (3000/100)  interferograms. 
The  bandwidth  is  increased  to  300  Hz. 

Interferogram  Averaging: 

The  averaging  of  interferograms  is  normally  accomplished 
by  sampling  at  a  rate  somewhat  greater  than  2fnax/  and  succes¬ 
sively  adding  the  values  to  a  series  of  partial  sums  stored 
in  a  computer  memory.  The  analog  to  digital  converter  must 
have  a  dynamic  range  of  at  least 

log 2  £lQ/N.  E.  P„  (kN/T)*5 

bits.  Since  k  scans  are  to  be  added  the  computer  word  length 
must  be  at  least 

log2  ^k5  Iq/N.E.P.  (N/T)*5] 
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bits,  although  only  the 


-°g. 


[ 


I  /N.E.P.  (N/T)' 


most  significant  bits  need  be  retained. 
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SECTION  IV 


DIFFERENTIAL  TECHNIQUES 

Differential  Interferograms : 

The  large  dynamic  range  typical  of  interferograms  re¬ 
sults  from  the  fact  that  I  is  given  by 

1 0  =  /  s<v>  dv 

AV 

where  S(v)  is  the  spectral  intensity  of  the  observed  source. 
I  can  be  greatly  reduced  by  causing  the  interferometer  to 
measure  the  difference  between  the  observed  source  and  a 
known  one.  This  procedure  is  facilitated  by  the  use  of  a 
dielectric  beam  splitter1  in  the  interferometer.  After  a 
discussion  of  the  theory,  several  examples  of  practical  in¬ 
terest  will  be  considered. 

Complementary  Outputs: 

If  the  interferometer  shown  in  figure  3  is  lossless, 
the  powers  of  the  reflected  and  the  transmitted  outputs  must 
sum  to  the  original  input  power.  Therefore,  any  modulation 
in,  say,  the  transmitted  light  caused  by  scanning  the 
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interferometer  must  be  accompanied  by  a  complementary  mod¬ 
ulation  in  the  reflected  light  so  that  energy  is  conserved. 
Figure  3  illustrates  this  situation  for  a  monochromatic 
input.  For  an  arbitrary  spectr,''  input  the  interf erograms 
.  ay  be  expressed  by 


VT> 


/ 


S(v)  cos  (2ttvt)  dv 


AV 


(8-a) 


IR(T) 


S  ( v)  cos  (2ttvt  +  tt)  dv 


(8-b) 


Dual  Inputs; 

The  interferometer  apertures  which  served  above  as 
outputs  may  each  simultaneously  act  as  an  input  to  the  inter¬ 
ferometer.  Suppose,  for  example,  that  the  interferometer  is 
illuminated  by  sources  S.^  and  S2  located  in  opposite  apertures. 
Since  the  interferometer  is  symmyetr ical  each  source  produces 
both  a  reflected  and  a  transmitted  interferogram  according 
to  (8).  The  outputs  each  consist  of  composite  interf erograms 
according  to 


=  J1T  + 


2R 


(9-a) 


I 


2 


(9-b) 
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REFLECTED 


Fij.  3  Co*ple*©nt*ry  ©stpats  of  lost  lets 

Miehelson  interf©roneter . 
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Due  to  the  symmetry  of  the  .interferometer  and  since  the  form 
of  (8)  does  not  depend  on  S,  we  can  write 


[1(V)  =  f  S1(V) 


/ 


Av 


v)  cos  2ttvt  dv  +  /  S2(v)  cjs(2ttvt  +  7 r)  dv 

Av 


j'  [Si^V^  ”  S2(v)j  COS  27rVT  dv* 


(10) 


Av 


The  second  output  I 2  is  -I  . 

From  (10)  we  see  that  the  dynamic  range  of  an  interfero- 
gram  can  be  arbitrarily  reduced  by  placing  an  appropriate 
known  source  S2  in  the  vicinity  of  the  detector.  This  source 
should  resemble  as  closely  as  possible  the  observed  source 
S^.  Our  experiment  would  therefore  consist  of  measuring  the 
difference  between  and  our  best  guess  of  S^.  The  dynamic 
range  requirements  are  reduced  according  to  our  a  priori 
information. 

The  use  of  equation  (1)  for  this  situation  is  conditional. 

If  (Sj  -  S2)  >  0  or  if  (S^  -  S2)  <  0  for  the  entire  spectrum, 

equation  (1)  is  applicable  and  S/m  indicates  the  mean 

s 

singal- to-noise  ratio  in  the  difference  spectrum  and  I  re- 

o 

lates  to  equation  (10).  The  final  spectrum  is  obtained  by 
adding  (or  subtracting)  S2  to  the  measured  difference  spectrum. 
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Equation  (1)  may  in  practice  be  applied  whenever  the  dif¬ 
ference  spectrum  is  predominently  of  one  sign.  When  this  ir 
not  the  case,  I  need  not  be  the  largest  portion  of  the 
inter ferogram  and  the  p'-ak  signal-to-noise  ratio  becomes  a 
function  of  the  form  of  and  and  is  therefore  difficult 
to  predict.  The  data  reduction  for  this  situation  is  compli¬ 
cated  by  the  fact  that  a  phase  sensitive  Fourier  transforma¬ 
tion  is  required  to  distinguish  the  negative  and  positive 
regions  in  the  difference  spectrum. 

Natural  Differential  Inter feroarams ; 

In  the  thermal  region  cf  the  infrared  the  detector  and 
its  environment  provide  a  second  input.  It  consists  of 
thermal  radiation  and  is  present  whenever  cooling  of  the 
instrument  is  not  performed  to  suppress  it.  The  instrument 
radiance  should  approx viia to  that  of  the  source  to  significant¬ 
ly  reduce  dynamic  range,,  For  convenience  where  two  sided 
inter ferograms  are  employe"5  -  phase  sensitive  Fourier  trans¬ 
formation  can  be  avoided  by  insuring  that  the  instrument 
spectral  radiance  is  always  somewhat  greater  than  that  of 
the  source. 

The  use  of  this  method  need  not  be  limited  to  the 
thermal  region.  In  spectral  regions  where  natural  thermal 
radiation  does  not  suffice,  the  second  input  can  consist  of 
a  hot  source  placed  in  the  neighborhood  of  the  detector. 
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SECTION  V 


SIGNAL  COMPRESSION 

Logarithmic  Amplification: 

The  most  plebian  approach  to  dynamic  range  problems  is 
logarithmic  signal  compression.  The  fact  that  i  iterferograms 
involve  both  positive  and  negative  signals  inflicts  only 
minor  variations  of  the  procedure.  For  example,  the  cube 
(or  fifth,  etc.)  root  of  the  signal  can  be  employed  in  place 
of  the  logarithm.  The  plot  of  these  curves  in  figure  4  shows 
that  the  desired  comprerssion  of  large  signals  is  retained, 
but  with  a  positive-negative  facility. 

Two  basic  problems  of  differing  severity  are  encountered. 
The  first  is  that  the  requisite  bandwidth  for  faithful  repro¬ 
duction  is  substantially  increased  unless  the  implementation 
is  carried  out  with  extraordinary  care.  The  effect  corresponds 
to  the  extra  harmonics  generated  by  noti-j  inear  distortion  in 
high  fidelity  (music)  systems.  Thus,  if  the  compression  is 
done  blithely  with  analog  circuitry  the  original  sample  fre¬ 
quency  will  be  inadequate.  More  exotic  techniques  such  as 
digital  manipulation  after  sampling  are  required  to  avoid 
the  problQm0 

The  second,  and  ultimately  more  severe,  problem  is  that 
the  errors  become  concentrated  in  regions  of  high  signal  of 
the  inter ferogram.  It  is  just  those  regions  which  provide 
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Fi«i.  4  Analytic  iifnal  eenpression. 


Fl§.  !S  Linear  eefnentei  eijnal  eenpretsion. 
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the  most  definitive  character  to  the  final  spectrum.  The 
space  frequency  distribution  of  the  noise  in  the  final 
spectrum  is  no  longer  white,  but  closely  corresponds  to  the 
space  frequency  distribution  of  the  spectrum.  The  presence 
of  such  non-uniform  noise  can  be  very  deceptive  for  the 
scientific  interpretation  of  the  spectrum.  Typical  character¬ 
istics  of  this  sort  will  be  evident  in  the  subsequent  section 
describing  experimental  (computer  simulated)  results. 

A  plausible  criterion  for  the  effectiveness  of  compres¬ 
sion  (or  other)  procedures  may  be  found  in  the  resultat 
distribution  function  of  the  interferogram.  One  should  strive 
to  utilize  each  level  of  the  transmission  channel  with  equal 
probability.  In  this  manner  each  word  conveys  a  maximum  of 
information. 

From  this  point  of  view,  the  problem  with  the  initial 
interferogram  is  that  it  tends  to  have  large  cr.etral  fringes 
but  then  spends  most  of  its  time  fluctuating  near  zero  level. 
The  higher  levels  of  the  transmission  channel  tend  to  be 
wasted  except  for  the  brief  cccurance  of  the  central  fringes. 

Gain  Changing; 

3 

Gain  changing  procedures,  such  as  advocated  by  Forman 
whereby  the  gain  is  decreased  when  the  signal  exceeds  a  pre¬ 
scribed  level,  are  simply  particular  cases  of  compression. 

In  this  case  we  obtain  a  characteristic  as  shown  in  figure  5 
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along  with  the  cubic  curve.  It  amounts  to  a  sectional  ap¬ 
proximation  of  the  smooth  procedure,  ant  engenders  the  same 
problems.  Torman's  results  clearly  show  the  non-white 
character  of  the  noise  on  is  resultant  spectra.,  although 
he  does  not  mention  the  effect. 

In  evaluating  this  procedure  it  must  be  kept  perfectly 
clear  that  one  bit  of  information  is  required  to  establish 
which  gain  pertains0  In  other  words,  the  spt.n  must  be  the 
total  vertical  span  of  the  linear  sections  and  not  just  the 
span  of  any  individual  section. 

Often,  one  can  avoid  this  "extra"  bit  with  a  modest 
amount  of  a  priori  information  concerning  the  spectrum.  Since 
large  amplitudes  occur  only  in  the  central  portion  of  most 
interferograms  it  is  possible  to  prescribe  the  lower  gain 
over  a  prespecified  central  portion  of  the  interferogram. 

The  a  priori  information  serves  to  ascertain  just  how  wide 
a  portion.  In  addition  to  eliminating  the  extra  bit,  this 
procedure  at  least  gives  consistently  non-white  noise  charac¬ 
teristics  which  are  signal  independent. 

We  have  also  considered  digital  data  transmission  systems 
involving  variable  word  lengths.  This  concept  led  to  such 
extensive  buffering  problems  however,  that  we  abandoned  it. 


Experimental  Simulations; 


A  series  of  computer  simulations  was  carried  out  on 
order  to  clearly  see  the  effects  of  dynamic  range  probxems. 

The  spectrum  shown  in  figure  6  was  adopted  and  its  cosine 
Fourier  transform  also  shown  in  figure  6  represents  the 
inter ferogram.  Next,  a  controlled  amount  of  random  noise  was 
added  to  the  interferogram,  and  the  result  retrans formed  into 
the  spectral  domain. 

The  effects  of  60,  50,  40,  and  30  dL  (voltage)  channels 
(Smax/ns  =  1000,  316,  100,  32.  respectively)  for  i nterferogram 
transmission  were  thus  simulated.  The  inter ferogram  from  the 
60  dB  channel  shown  in  figure  7  appears  to  the  eye  as  a  faith¬ 
ful  rendition  of  the  original.  The  resulting  spectrum  shown 
in  figure  7  is  conspicuously  degraded  however.  That  comparison 
serve.s  to  illustrate  the  severity  of  the  dynamic  range  prob¬ 
lem.  Subsequent  figures  for  the  50,  40,  and  30  dB  channels 
show  that  when  noise  becomes  conspicuous  on  the  interferogram, 
the  spectrum  becomes  barely  recognizable. 

Turning  to  the  compression  procedures,  figure  11  shows 
the  cube  root  compressed  interferogram.  The  faint  trace  is 
the  original  interferogram  superimposed.  Noise  simulating 
a  40  dB  channel  was  added  to  the  compressed  interferogram, 
the  interferogram  expanded  and  the  spectrum  shown  in  figure 
11  resulted.  On  the  face  of  it,  this  spectrum  seems  comparable 
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Spectrum  No  Noise 


Pif.  6  Interferogran,  Spectrun 
adopted  for  sinulationa. 
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Fifl*  ?  Interf erof raw,  Spectrum 
ewployinj  60  <B  channel. 


Fiji,  8  Intsrf erogram,  Spectrum 
employing  50  channel. 
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to  a  55  dB  channel  without  compression,  and  so  we  estimate  a 
15  dB  decrease  in  dynamic  range  requirements. 

That  conclusion  is  very  deceptive.  Note  in  particular 
the  wandering  zero  level.  This  w'ould  not  be  at  all  apparent 
without  the  noiseless  spectrum  for  comparison.  The  wandering 
level  is  due  to  the  non-white  character  of  the  noise,  and  re¬ 
flects  the  ,'upplementary  low  space  frequency  noise. 

Already  in  the  interferogram  an  incipient  fault  of  the 
procedure  is  becoming  evident.  This  fault  is  that  the  signal 
is  displaying  a  marked  avoidance  of  zero  level.  The  effect 
would  be  even  more  pronounced  for  fifth  root  compress icn. 

The  consequence  of  the  effect  is  that  the  transmission  chan¬ 
nel  is  still  used  inefficiently.  The  signal  distribution 
function  greatly  upsets  uniform  probability  of  the  available 
channel  levels. 

The  simulation  of  a  sectional  gain  change  is  shown  in 
figure  12.  The  gain  of  the  outer  fringes  was  increased  by  a 
factor  of  32  (30  dB) ,  then  the  noise  corresponding  to  a  30  dB 
channel  added.  The  interferogram  shown  in  figure  12  was  re¬ 
cons  tr  cted  from  the  result.  This  looks  like  a  30  dB  channel 
for  the  central  fringes  and  the  equivalent  of  a  60  dB  chann  1 
for  the  outer  fringes.  The  compressed  interferogram  could 
have  been  genuinely  transmitted  on  a  30  dB  channel. 
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2  Channel  Interferogram  Gain 
Ratio  32  30  dB  Channel 


Spectrum  2  30  dB  Channels 
8  to  1  Gain  Ratio 


Fig.  12  Interferogram,  Spectrum 
employing  factional  gain  change. 
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Again,  the  noise  seems  to  correspond  to  a  15  or  even 
20  dB  improvement,  but  also  again,  the  wandering  zero  level, 
or  low  space  frequency  noise  is  evident. 
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SECTION  VI 


CHIRPING  TECHNIQUES 


Chirping: 

Horizontal  expansion  of  the  central  fringes  offers  an 
alternative  to  the  vertical  compression  techniques.  Horizon¬ 
tal  expansion  comes  from  making  the  position  of  the  central  fringes 
frequency  dependent.  The  resultant  fringes  resemb]  i  chirp 
interferogram  is  not  inappropriate. 

Although  the  immediate  present  objective  of  this  study 
of  the  chirp  procedure  is  the  reduction  of  dynamic  range,  a 
more  far  reaching  consequence  is  the  additional  insight  i 
provide  5  concerning  the  nature  of  inter ferograms. 

Shape  of  the  Central  Fringes: 

A  cosinusoidal  form  is  usually  adopted  to  describe  the 
fringes  of  a  Michelson  interferometer.  Thus,  the  fringes 
are  symmetric  about  zero  path  difference  and  the  fringes  from 
all  wavelengths  are  all  in  phase,  constructively  interfering 
to  provide  a  large  amplitude  fringe.  Phase  shifts  upset 
this  description  with  varying  degrees  of  severity.  A  few 
simple  examples  are  in  order. 

First,  the  introduction  of  a  phase  shift  independent 
of  frequency  destroys  the  symmetry  of  the  interferogram,  while 
retaining  the  envelope  of  the  interferogram.  For  example, 
a  90°  phase  shift  converts  perfect  symmetry  to  perfect: 
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antisymmetry.  It  is  important  to  realize  that  any  angle  may 
be  invoked  so  that  the  resultant  inter ferogram  will  generally 
be  non- symmetric,  iven  in  this  situation  all  fringe  frequen¬ 
cies  have  the  same  phase  at  tha  center,  and  therefore  add 
constructively.  It  is  very  important  to  be  aware  of  the  ex¬ 
istence  of  such  a  non-symmetric  but  unchirped  interferogram. 
In  normal  Fourier  spec trome ter y  this  is  a  frequent  experimen¬ 
tal  occurance  and  should  not  be  disconcerting. 

For  a  second  example,  the  introduction  of  a  phase  shift 
which  is  proportional  to  fringe  frequency  serves  only  to  dis¬ 
place  the  entire  interferogram.  The  form  of  the  central 
fringes  remains  unaltered,  and  the  simple  shift  of  our  frame 
of  reference  restores  the  original  symmetry. 

Any  more  complicated  dependence  of  phase  on  fringe  fre¬ 
quency  can  be  thought  of  in  terms  of  linear  segments.  Each 
linear  segment  would  require  a  different  shift  in  the  frame 
of  reference  according  to  the  slope  of  that  segment.  In 
this  fashion  the  central  fringe  position  may  be  distributed 
as  a  function  of  fringe  frequency  over  a  large  retardation. 
The  total  "power"  in  the  interferogram  remains  unaltered  and 
so  the  maximum  fringe  amplitude  diminishes.  Roughly  speaking 
if  we  quadruple  the  half  width  of  the  central  finges,  the 
peak  amplitude  is  halved  (square  root  relationship). 
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The  notion  of  stationary  phase,  rather  than  zero  path 
difference,  now  serves  to  specify  the  central  finges.  AS 
long  as  the  phase  is  a  smooth  function  of  finge  frequency 
then  an  appropriate  shift  of  the  reference  frame  can  nul 
the  gradient  at  any  one  frequency.  For  this  shifted  position 
the  phase  is  independent  of  frequency  for  a  small  interval, 
and  the  center  of  the  reference  frame  suffices  to  dfine  the 
central  finge  position  for  that  frequency  interval. 

Any  more  detailed  analysis  would  quickly  reveal  that  this 
description  of  fringes  amounts  to  group  velocity  propagation. 
The  phase  shifts  are  most  easily  introduced  by  the  insertion 
of  refractive  material  in  to  one  beam  of  the  interferometer. 

The  phase  shifts  result  from  dispersion  in  the  index  of  re¬ 
fraction  (phase  velocity)  of  the  material.  The  position 
shifts  of  the  '.entral  finges  can  be  descr  bed  in  terms  of 
group  velocity  index  of  refraction  of  the  material. 

The  actual  details  are  not  particularly  relevant  in 
this  report.  Suffice  it  to  say  that  two  different  materials, 
one  in  each  beam,  must  be  used  if  interferometer  throughput 
is  to  maintained.  Detailed  formulae  for  the  appropriate  thick¬ 
ness  may  be  found  in  reference  1. 

Experimental  Simulation: 

The  same  test  spectrum  served  for  a  computer  simulation 
of  chirping.  For  a  small  touch  of  realism,  the  spectrum  was 


33 


considered  to  range  from  about  4.5  to  0.8  micron,  with 
the  dispersive  effects  of  5  mm  of  IRTRAN-5  (MgO)  pi  viding 
the  chirping.  Figure  13  shows  the  chirped  and  unchirped 
inter ferogr am.  The  central  finge  position  in  this  case  is 
not  a  monotonic  function  of  frequency  but  bends  back.  Thus, 
low  frequencies  are  on  the  right,  intermediate  at  the  left, 
and  the  high  frequencies  return  a  little  toward  the  right. 

The  reduction  of  dynamic  range,  found  from  the  ratio 
of  the  maxima  in  the  chirped  and  unchirped  interferograms 
is  about  13  dB.  The  spectrum  calculated  from  the  chirped 
inter ferogram  is  shown  in  figure  13.  In  t1  is  case,  the 
amplitude  and  phase  Fourier  transform  served  in  place  of  the 
previous], y  employed  cosine  transforms.  There  was  clearly 
some  malfunctioning  of  the  program  or  plotter,  but  the  re¬ 
sults  looked  sufficiently  encouraging  to  do  some  real  experi¬ 
ments. 

In  order  to  appreciate  the  remaining  “xperimenta.l  results, 
a  discussion  of  the  possible  data  reduction  procedures  is 
essential 

On  The  Reduction  of  Extremely  Chirped  Interferograms: 

An  extremely  chirped  incerf erogram  is  one  in  which  the 
distribution  of  central  finge  positions  as  a  function  of  color 
spans  almost  almost  the  entire  length  of  the  inter ferogram. 

It  is  even  permitted  to  fold  back  so  that  both  short  and  long 
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Fig.  13  Interfer*gra«,  Spectra* 
enploying  chirping. 


wavelengths  have  their  central  fringes  near  one  end  of  the 
inter ferograin  while  the  middle  wavelengths  have  their  cen¬ 
tral  fringes  near  the  other  end.  We  will  subsequently  see 
that  if  the  chirping  is  restricted  to  the  central  two  thirds 
of  the  interferogram  the  analysis  remains  convenient.  It 
will  be  shown  in  the  appendix  that  extreme  chirping  is  still 
sufficiently  mild  to  satisfy  the  conditions  of  the  general 
relation^  Thus,  as  long  as  a  phase  curve  is  available  (this 
must  be  obtained  as  the  phase  spectrum  of  a  smooth  black 
body  spectrum) ,  phase  correction  is  legitimate. 

The  clumsy  aspect  of  straightforward  phase-corrected 

Fourier  analysis  is  that  the  instrumental  profile  is  misleading 

for  those  spectral  regions  whose  central  finges  lie  away  from 

the  center  of  the  interferogram.  The  clumsiness  is  actually 

4/ 

more  apparant  via  the  spectral  transfer  function  — '  (the 
cosine  Fourier  transform  of  the  instrument  profile).  Accor¬ 
ding  to  the  general  relation  the  spectral  transfer  function 
is  just  the  symmetric  part  of  the  truncation  function  about 
the  central  fringe. 

Figure  14  serves  as  a  guide  tc  figures  15  through  18,  and 
shows  the  location  of  the  central  fringes  (vertical  slash) 
with  respect  to  the  interferogram  range  (boxcar)  for  a  through 
_i.  Figure  15  shows  the  spectral  transfer  functions  appli¬ 
cable  to  straightforward  calculation  (phase-corrected  Fourier 
transformation)  for  the  range  of  central  fringe  locations 
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Central  fringe  peaition  relative  te  truncation  boxcar  for  various  offeete 
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offsets  with  trian^slar  apsdisation* 
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Speetral  transfer  functions  for  various  offsets 


£  through  i.  .  When  the  central  fringes  are  not  centrally  lo¬ 
cated  in  the  truncating  function,  and  ungainly  steo  occurs 
in  the  spectral  transfer  function.  A  step  discontinuity 
leads  to  sidelobes  on  the  instrument  profile,  which  decay  as 

1/Af,  whereas  a  slope  discontinuity  leads  to  a  decay  rate  of 
2 

1/Af  .  When  the  step  appears  as  in  figure  15  _f  the  contrast 
of  unresolved  absorption  lines  is  diminished  by  about  1.6 
from  that  which  one  would  anticipate  based  on  instrumental 
resolution. 

J.t  would  be  possible  to  resort  to  two  dimensional  apodiza- 
tion  (apodization  as  a  function  of  wavelength  as  well  as  re¬ 
tardation)  ,  but  this  precludes  the  use  of  the  efficient 

6/ 

Cooley-Tukey  transform  program  Much  less  drastic  measures 

are  definitely  called  for,  and  will  be  perfectly  suitable. 

If  triangular  apodization  is  applied  to  the  inlerfero- 
gram,  the  resulting  set  a  spectral  transfer  functions  is  shown 
in  figure  16.  These  functions.-  at  least  _a  through  are 
obviously  superior  to  those  of  figure  15.  The  only  defect  is 
easily  removed  by  equalizing  the  spectrum.  A  given  wavelength 
has  its  central  fringes  a  certain  amount  off-center  (as  deter¬ 
mined  from  the  black  body  re Cerence  phase  spectrum)  and  so  must 
be  multiplied  by  a  factor  to  normalize  that  spectral  transfer 
function.  The  anodization  remedies  the  contrast;  the  equal¬ 
ization  restores  the  instrumental  wavelength  response.  Two 
other  initial  apodizations  are  shown  in  figures  17  and  18. 
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The  first  is  a  half  cosine  wave,  the  second  a  trapezoid, 
appearing  as  in  17  and  18  respectively.  Clearly  these  apo- 
dizations  are  suitable  out  to  offsets  between  f  and  £.  Thus, 
if  the  chirping  ' s  restricted  to  roughly  the  central  two 
thirls  of  the  interferogram,  this  apodization-equal ization 
procedure  is  very  suitable.  For  larger  offsets,  the  hump 
appearing  in  the  spectral  transfer  function  will  become  un¬ 
acceptable,  and  the  last  offset  _i  (to  the  edge)  is  impossible 
to  equalize. 

Offset  in  the  range  £  to  h  can  be  mac  e  available  by  using 
a  second  pass  of  the  Cooley-Tukey  program  using  a  different 
initial  apodization  and  then  linearly  combining  the  two  spec¬ 
tra  with  weighing  fact  s  dependent  on  wavelength.  Cursory 
investigation  of  this  two  (or  more)  pass  composite  procedure 
indicates  that  the  rewards  are  probably  not  worth  the  trouble. 

We  may  no  v  recapitulate  the  recommended  computing  pro¬ 
cedure.  A  reference  interf  •rogram  is  obtained  of  a  continuum 
source  such  as  a  black  body.  This  interferogram  should  then 
be  triangularly  apodized  (apodization  is  desirable  but  not 
essential)  and  the  amplitude  and  phase  spectra  computed. 

These  will  appear  as  in  figure  18A,  phase  always  being  referred 
to  the  central  sample.  With  the  understanding  that  the  phase 
spectrum  is  modulo  2tt,  the  step  discontinuities  are  removed 
by  adding  integral  multiples  of  2tt  (actual  computation  will 
probably  be  in  cycles  rather  than  radian;  1  cycle  -  2ir  radians). 
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A  least  squares  polynomial  is  then  fit  to  the  phase  spectrum. 
This  polynomial  is  then  differentiated.  This  derivative  is 
a  measure  of  the  offset  cf  the  central  finges  versus  fre¬ 
quency.  Its  scale  factor  is  one  cycle  per  sample  frequency 
per  sample  point. 

The  unknown  interf erogram  is  next  apodized  (triangular, 
cosine  or  trapezoid)  and  Fourier  transformed.  From  its 
phase  spectrum,  we  subtract  the  curve  fit  of  the  reference 
phases.  The  result  should  be  an  approximately  straight  line, 
almost  horizontal  and  almost  intercepting  zero.  J.f  it  is 
not  accurately  horizontal  and/or  does  not  intercept  zero, 
it  should  be  modified  to  do  so  by  adding  a  constant  and  a 
gradient.  Phase  correction  is  now  carried  out  by  multiplying 
the  amplitude  spectrum  by  multiplying  the  ampltidue  spectrum 
by  the  cosine  of  the  modified  phase  spectrum. 

Equalization  is  n  xt  in  line.  The  reciprocal  of  the  equal¬ 
ization  factor  is  simply  the  apodization  factor  at  the  offset 
determined  two  paragraphs  ago.  If  any  appreciable  modification 
was  applied  to  the  phase  spectrum  by  adding  a  linear  term, 
the  corresponding  offset  modification  (a  constant)  should  pre¬ 
cede  ascertaining  the  apodization  (1/equalization)  factor. 

The  phase  corrected  spectrum  is  multiplied  by  the  equalization 
fac tor. 

This  final  spectrum  is  that  of  the  source,  with  spectral 
transfer  functions  as  indicated  by  the  figures.  In  other  words. 
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the  instrument  profile  changes  somewhat  over  the  spectrum  to 
correspond  to  the  spectral  transfer  function  at  the  offset 
corresponding  to  each  spectral  frequency. 

Txperimental  Results: 

The  real  experiments  involved  an  interferometer  chirped 
with  5  mm  of  IRTRAN-5  vs.  9  mm  of  CeiF^  in  the  other  beam. 

The  source  was  an  acetylene  flame,  and  the  wavelength  range 
was  mostly  between  2  and  2.5  microns.  Half  sample  frequency, 
the  rightmost  point  on  the  spectral  traces,  was  at  5870  cm  1. 

The  position  of  the  central  fringes  as  a  function  of 
wavelength  is  shown  in  Figure  19.  This  position  was  determined 
experimentally  as  the  slope  of  a  black  body  phase  spectrum. 

The  factor  which  relates  position  to  slope  is  one  sample  point 
for  each  cycle  per  sample  frequency. 

For  this  extreme  degree  of  chirping,  the  questions  of 
validity  of  data  reduction  procedures  are  paramount.  A 
chirped  acetylene  interferogram  in  shown  in  f.igure  19A.  The 
amplitude  and  sphase  spectra  of  a  similar  interferogram  is 
shown  in  figure  20.  It  is  interesting  to  note  the  resemblence 
of  the  envelope  of  the  interferogram  to  the  amplitude  spectrum. 
While  this  similarity  is  not  accidental,  it  is  of  little  use 
as  a  spectrum.  For  example,  monochromatic  light  still  pre- 
snets  uniform  fringes  over  the  entire  interferogram. 
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Several  procedures  were  used  in  the  data  reduction. 

Simple  amplitude  and  phase  spectra  served  for  figure  20. 

Phase  correction  led  to  figure  21.  Note  that  phase  correc¬ 
tion  he  -  a  hardly  noticeable  effect  of  the  spectrum. 

The  procedure  of  trapezoidal  apodization  followed  by 
equalization  led  to  figure  22.  Note  that  the  series  of  lines 
near  3400  cm  ^  have  vanished. 

Finally,  the  transform  was  carried  out  with  two-dimensional 
apodization  (fig.  23).  This  required  an  extremely  lengthy  calcu¬ 
lation  since  the  Cooley-Tukey^f actor ing  algorithm  was  not  applicable. 
Note  that  the  contrast  of  that  series  of  bands  is  doubled. 

The  band  contrast  is  readily  understood  in  terms  of  spec¬ 
tral  transfer  function.  A  uniform  band  series  like  that  gives 
a  blip  in  the  int.erf  erogram,  in  this  case  at  the  leftmost 
edge  of  figure  19.  The  central  fringes  for  this  spectral 
integral  lie  in  the  righthad  major  hump.  Thus,  the  corres¬ 
ponding  righthand  blip  due  to  the  band  lies  outside  the 
inter f erogram.  In  figure  24  we  see  the  three  pertinent 
spectral  transfer  functions  and  the  location  of  the  blip  is 
indicated. 

No  apodization  leads  to  a  contrast  transfer  of  \  at  the 
space  frequency  of  the  blip.  Simple  trapezoidal  apodization 
leads  to  almost  zero  contrast  transfer.  Full  two-dimensional 
apodizing  leads  to  unit  contrast  transfer.  The  experiments 


47 


ACETYLENE 

ONE  SCAN  INTERFEROGRAM 
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Fig.  19A  Chirped  interferogran 
of  acetylene  flane  source. 
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20  Amplitude,  Phase  spectra 
no  apod izat ion. 
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Pig.  21  Phase  corrected  spectrum. 
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Fig.  22  Spectrum  from 
trapezoidal  apodization,  equalized. 
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.  23  Spectrun  fron 
i«ensional  apodization 
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Pip.  24  Spectral  transfer  fanettona  for  bania 
near  3400  cn”1. 


54 


agree  perfectly  with  expectations,  and  serve  to  highlight  the 
attributes  of  the  various  procedures. 

Phase  Correction: 

1  7  8/ 

Normally,  phase  correction  procedures  — 1  J  f  obtain  a 
reference  phase  spectrum  from  the  narrow  group  of  central 
fringes.  Peripheral  fringes  are  excluded  since  the  reference 
phase  spectrum  is  known  to  be  smooth  (without  high  resolution 
detail)  and  so  the  peripheral  fringes  would  only  contribute 
noi.se. 

With  extremexy  chirped  interferograms  there  is  no  narrow 
group  of  ventral  fringes  and  we  are  obliged  to  resort  to  another 
approach.  This  approach  is  to  calculate  the  phase  spectrum 
from  a  black  body  source,  and  use  this  as  a  reference  phase 
spectrum  after  suitable  fitting.  Suitable  fitting  involves 
adjusting  only  two  parameters;  the  constant  and  linear  ferns 
of  the  reference  phase  spectrum.  The  higher  order  terms  are 
reproducible. 

It  turns  out  to  be  easy  enough  to  accomplish  a  visual 
fit;  automating  the  process  is  more  difficult.  The  problem 
is  akin  to  character  recognition.  None  of  our  attempted 
methods  functioned  with  complete  consistency  and  many  did  not 
function  at  all,  probably  due  to  residual  program  bugs.  Cer¬ 
tain  aspects  are  nevertheless  worthy  of  elaboration. 
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Two  problems  and  the  combination  thereof  prove  paramount.. 
One  is  360°  discontinuities  encountered  since  the  phase  is 
mathematically  modulo  2ir ,  but  physically  continuous  over 
many  cycles.  Such  discontinuities  must  be  eliminated  if  the 
fitting  of  the  reference  phase  is  to  be  at  all  meaningful. 

The  second  problem  is  the  wild  fluctuations  of  phase  wherever 
the  signal- to-noise  ratio  is  loss  than  unity.  These  fluctua¬ 
tions  have  no  bearing  on  the  fit  and  must  be  properly  dis¬ 
criminated  against,, 

As  far  as  the  elimination  of  the  discontinuities  is 
concerned,  two  methods  exist.  The  first  is  to  look  for  any 
steps  aiong  the  phase  spectrum  greaterri  than  w  in  magnitude. 
Appropriately  adding  or  subtracting  2w  thereafter  eliminates 
the  steps. 

The  discontinuities  may  be  also  eliminated  by  differen- 

9/ 

tiating  and  integrating  the  phase  spectrum  .  Since  the 
analysis  is  discrete  finite  differences  are  more  legitimate. 


= 


k+1 


-  <*V 


where  <I>^  is  the  phase  at  frequency  k.  C  and  S  are  cosine  and 
sine  terms  of  the  Fourier  transform. 


tanA$k 


C.  S,  -  S,  CL  , 
k  k+1  k  k+1 

C,  C,  .  +  S.  S,  , 
k  k+1  k  k+1 


(ID 
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It  is  interesting  that  when  Atf^  remains  small  then  tanA<t>^ 

«  A<J^,  and  the  approximate  phase  spectrum  may  be  genera¬ 
ted  as  successive  partial  sums  of  the  right  hand  side. 

A  very  useful  mean  slope,  to  determine  the  center,  of 
an  unchirped  inter ferogran  may  be  readily  obtained  from 


equation  (11).  Note  that  if  and  S^.  &  then 

2 

the  denominator  of  (11)  is  A  .  With  this  we  can  very  easily 


create  an  A  weighted  average  of  tanA$  as 


tanA<£>  = 


k  (CkSk-?-l  ~  SkCk+l) 

£  <ckck.i  +  vw 


from  which  we  can  readily  determine  the  appropriate  gradient 


A  $  =  arctan  (tanAO)  .  (13) 

This  procedure  is  very  simple  and  elegant  v  .id  certainly 
ought  to  work,  in  spite  of  the  fact  that  we  did  not  succeed 
with  it  in  practice. 

The  procedure  should  need  only  modest  variation  for 
chirped  inter ferograms ,  which  would  effectually  subtract  the 
higher  (than  first)  order  terms  of  the  reference  phase  before 
averaging. 

The  only  pru^adure  which  has  so  far  worked  in  practice 
is  abritratily  selecting  those  segments  of  the  spectrum  whose 
amplitude  everywhere  exceeds  10%  of  the  maximum  amplitude. 
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360°  discontinuities  were  eliminated  by  the  first  process. 

A  least  squares  line  was  then  fit  to  the  difference  between 
that  segment  of  the  phase  spectrum  and  the  reference  phase 
spectrum.  Differences  from  this  line  were  finally  used  as 
the  angle  deviations  for  phase  correction. 

Describing  the  Reference  Phase; 

In  order  to  provide  a  smooth  reference  phase  spectrum 
for  subsequent  computations ,  the  phase  spectrum  of  a  black 
body  source  is  obtained.  The  360°  discontinuities  are  elim¬ 
inated  and  within  the  frequency  range  of  good  signal-tc-  r~ise, 
a  least  squares  polynomiax  is  fit  to  the  curr 

In  practice,  with  extreme  ch-rping,  a  considerable  number 
of  terms  are  required  for  a  good  fit;  too  mai.y.  Considering 
the  chirping  range  and  required  accuracy  fc  situation  is 
not  surprizing.  For  the  interferometer  used  in  the  examples 
of  the  section  entitled  Experimental  Results,  the  phase  spec¬ 
trum  spanned  about  800  cycles.  Considering  that  the  reference 
phase  should  be  accurate  to  better  than  l/40th  cycle,  the 
polynomial  need  be  accurate  to  better  than  one  part  in  30,000. 

Polynomial  fits  of  even  ninth  order  did  not  give  suf¬ 
ficient  accuracy.  The  very  best  obtainable  accuracy  was 
about  15°.  The  error  was  clearly  systematic  and  furthermore 
higher  order  fits  proved  exceedingly  slow  in  reducing  the 
systematic  errors.  Thus,  ordinary  Taylor  series  are  clearly 
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inappropriate  for  this  situation. 

In  particular,  hyperbolic  character  is  bestowed  on  the 
phase  spectrum  by  the  dispersion  (theoretical)  of  the  chirp¬ 
ing  material.  Fitting  a  hyperbola  with  a  power  series  can  be 
discouraging. 

The  remedy  proved  to  be  to  start  the  Taylor  series  with 
negative  exponents.  In  this  fashion  the  error  was  reduced  to 
less  than  2°  using  the  -2nd  to  +5th  powers  of  the  frequency. 
That  accuracy  is  quite  sufficient. 
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SECTION  VII 


FINAL  COMMENT 


A  quotation  from  Fellgett  is  suitable  to  conclude 
this  report,  "A  remark  on  noise:  Dynamic  range  becomes 
a  problem  for  high,  not  for  low  S/N  ratio. " 
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APPENDIX 


Lat  us  adopt  N  input  samples  which  give  rise  tc  N  output 

frequencies.  The  proportionality  factor  in  the  phase  slope 

versus  offset  relation  is  one  cycle  per  sample  frequency 

par  sample  point.  Therefore,  if  the  phase  slope  is  zero  at 

one  end  of  the  interferogram,  it  is  unity  at  the  other  end. 

At  this  far  end  the  phase  changes  by  one  cycle  per  output 

point.  Distributing  this  slope  change  uniformly  over  the 

•  • 

ouput  points  gives  <J>  =  1/N. 

Typically,  the  sidelobe  structure  of  the  instrument 
profile  is  thoroughly  depicted  by  less  than  20  output  points. 
Thus,  if  N  =  1000,  the  maximum  phase  e.  ror  involved  is 
\ (20)/1000  =  1/100  wave.  Thus,  for  each  instrument  profile 
(both  its  symmetric  and  antisymmetric  components  as  used 
in  the  general  relation)  the  interferogram  satisfies  the 
stationary  phase  condition  within  roughly  1/100  wave. 
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